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TRAVELING GRAVITY WATER WAVES WITH CRITICAL 

LAYERS 


AILO AASEN AND KRISTOFFER VARHOLM 


Abstract. We establish the existence of small-amplitude uni- and 
bimodal steady periodic gravity waves with an affine vorticity distribution, 
using a bifurcation argument that differs slightly from earlier theory. The 
solutions describe waves with critical layers and an arbitrary number 
of crests and troughs in each minimal period. An important part of 
the analysis is a fairly complete description of the local geometry of the 
so-called kernel equation, and of the small-amplitude solutions. Finally, 
we investigate the asymptotic behavior of the bifurcating solutions. 


1. Introduction 

Up until fairly recently, most authors working with steady water waves 
have made the assumption that the vorticity 

( 1 . 1 ) uj :=Vx- Uy 

of the velocity held (n, v) vanishes identically. Such waves are known as 
irrotational, as opposed to rotational waves where oj is allowed to be nonzero. 
Rotational waves can exhibit more exotic behavior than irrotational ones, in¬ 
cluding interior stagnation points and critical layers of closed streamlines [10]. 
Stagnation points correspond to huid particles that are stationary with re¬ 
spect to the wave, and for irrotational hows this can only occur at a sharp 
crest [34]. 

Irrotational waves are mathematically simpler to work with than rotational 
ones, due to the existence of the velocity potential. The velocity potential 
is the harmonic conjugate of the stream function, thus enabling the use of 
tools such as complex analysis, which are typically not available with nonzero 
vorticity. The survey [31] treats the theory of Stokes waves—an important 
class of irrotational waves—and the results on the so-called Stokes conjecture 
for such waves. This conjecture was not fully settled until the appearance of 
the paper [27]. 

Although rotational waves were considered intractable for mathematical 
analysis, they have long been important in more applied helds because 
rotational waves are not uncommon in nature: There are many physical 
effects that can induce rotation in waves, such as wind and thermal or salinity 
gradients [25], and rotational waves are also important in wave-current 
interactions [30]. 
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The first, and still the only known, explicit example of a nontrivial traveling 
gravity water wave solution to the Euler equations was given in [14] (see also 
[2] for a more modern treatment) and is rotational; a fact which was only 
later pointed out by Stokes. Much later came the first existence result for 
small-amplitude waves with general vorticity distributions [9]. It was not, 
however, before the pioneering article [4] that large-amplitude waves were 
constructed, using an extension of the global bifurcation theory of Rabinowitz 
[28, 16], leading to renewed interest in rotational waves. A corresponding 
result on deep water, where the lack of compactness is an obstacle, was 
established in [18]. 

Due to the methods used, neither the waves in [9] nor those in [4] exhibit 
stagnation. The first waves with a critical layer were constructed in [35], 
having constant vorticity. A different approach was used in [7], allowing for 
wave profiles with overhang (for which existence is still an open question, 
with some numerical evidence in the affirmative [32]). The method of proof 
for the existence of nontrivial rotational waves is typically bifurcation from 
parallel flows with a prescribed vorticity distribution. Such parallel flows are 
described in great detail in [22]. 

Other authors have looked at waves with density stratihcation [13, 17, 36], 
waves with compactly supported vorticity [29, 33], waves with discontinuous 
vorticity [5], and waves with a general vorticity distribution and stagnation 
[23]. An upcoming result also establishes the existence of large-amplitude 
gravity water waves with a critical layer [6]. This was done in the presence 
of capillary effects in [24], using an entirely different formulation. 

Of particular interest to us are [11, 12], which cover small-amplitude waves 
with an affine vorticity distribution. This is the natural step up from the 
constant vorticity considered in [35], and the resulting waves can have an 
arbitrary number of critical layers [10]. 

In this paper, which builds upon [1], we consider the same setting as 
in [11]. Small-amplitude solutions with an affine vorticity distribution are 
found by bifurcating from trivial solutions that depend naturally on three 
parameters. By using other choices for the bifurcation parameters in our 
argument, we obtain solution curves and sheets that, in general, do not 
coincide with those found in [11]. We are led to examine the asymptotic 
behavior of the bifurcating solutions; in particular for carefully chosen special 
cases. A complicating factor for our choice of bifurcation parameters is that 
they require an additional condition on the parameters. This condition can 
be interpreted as a nondegeneracy condition for the equation governing the 
dimension of the linearized problem. 

Another novel aspect of this work is a fairly complete description of the 
local geometry of the kernel of the operator appearing in the linearization. 
This is used to describe the geometry of the solution set near any trivial 
solution where the linear problem is one-dimensional, and for a class of 
trivial solutions with a two-dimensional linearized problem. We also show, by 
explicit construction, that the dimension of the linear problem can become 
arbitrary large for certain wavenumbers. This opens up the possibility for 
waves with arbitrarily many modes. Finally, we prove a regularity result, 
showing that the solutions we hnd are real analytic. 
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The outline of the paper is as follows: In Section 2 we formulate the problem 
and describe the setting in which we will work. Next, Section 3 focuses on the 
kernel of the linearized operator. Section 4 contains the bifurcation result for 
a one-dimensional kernel and gives the properties of the resulting bifurcation 
curves, while the final section. Section 5, covers two-dimensional bifurcation. 
Some useful derivatives are listed in Appendix A. 

2. The governing equations 

We consider pure gravity waves. The fluid motion is assumed to be 
incompressible and two-dimensional, with the coordinate system oriented so 
that the x- and y-axes are horizontal and vertical, respectively. The fluid 
domain is bounded below by a flat bottom, and above by a free surface. 
Within this setting, our aim is to construct solutions of the steady water-wave 
problem; that is, to hnd a surface prohle r] and a velocity held {u, v), dehned 
in the huid domain 

Qr/ ■= {(x, y) G : 0 < y < d + r/(x)j, 
where d is the depth of the undisturbed huid, satisfying the Euler equations 
(2.1a) Ux + Vy = 0, 

(2.1b) {u-c)ux + vuy =-px, 

(2.1c) (u - c)vx + vvy = -py - g 

in The surface prohle is assumed to satisfy r] > —d, so that the bottom 
is not exposed to air. In (2.1b) and (2.1c) the quantity p is the pressure, g 
is the gravitational acceleration, and c is the constant velocity at which the 
wave travels. 

In addition to the equations in (2.1), we impose the boundary conditions 


(2.2a) 

u = 0 

at y = 0. 

(2.2b) 

H 

T? 

1 

II 

at y = d 

(2.2c) 

p = 0 


The hrst two boundary conditions are known as kinematic boundary condi¬ 
tions, and state that there is no hux through the surface or bottom. The 
dynamic boundary condition in (2.2c) ensures that there is no jump in 
pressure across the free surface. 

We will be searching for periodic waves only, and so we introduce the 
wavenumber k > 0, and stipulate that all functions above be 27r/K-periodic 
in the horizontal variable. 

2.1. Stream function formulation. We now reformulate the water wave 
problem (2.1)-(2.2) in terms of a potential i/:, called the relative stream 
function. From incompressibility (2.1a), together with Qy being simply 
connected, we know that there exists a function f;: Qy ^R satisfying 

xfx = -V, ify = u - C. 

This function is uniquely determined by {u,v), up to a constant. 

The kinematic boundary condition (2.2a) is equivalent to '0a; = 0 at y = 0, 
and so is constant on the bottom. Similarly, we can use (2.2b) to deduce 
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that ijj is constant also on the surface. Next, (2.1b) and (2.1c) can be used 
to show that 

{V’, AV’} = 0, 

where {•, •} is the Poisson bracket defined by 

{fid} fydx ~ fx9y 

Furthermore, by also using the boundary conditions in (2.2), one can infer 
that the surface Bernoulli equation 

+ OT = Q on y = d + ri{x) 

holds for some Q G M. 

In terms of the stream function, the vorticity is given by 

oj = —A'0, 

which follows directly from its dehnition in (1.1). Observe also that is 
27r/K-periodic in the horizontal variable. To see this, note that {x,y) i—)• 
V’(x + 2TrfK, y) is also a stream function, taking the same constant values as 
■0 on the boundary. By uniqueness, they must be identical. 

The motivation for introducing the stream function is that, for a pre¬ 
scribed vorticity, the preceding equations are in fact equivalent to the steady 
water-wave problem. A precise statement, taken from [10], can be found in 
Proposition 2.1 below. We will use a subscript k, to denote 27r/K-periodicity 
in the horizontal variable. 


Proposition 2.1 (Stream function [10]). For y G C'^(M), u,v G (7^(11^) and 
a prescribed vorticity oj G (7^(11^), the steady water-wave problem (2.1)-(2.2) 
is equivalent to the stream function formulation 


(2.3a) 


A0 = —uj 
{0, A0} = 0 

0 = mo 
xjj = mi 

+ gy = Q 


in Qri, 

at y = 0, 

at y = d-\- r]{x), 


for 0 G C^{Qrj) and constants mo, mi and Q. 


2.2. The vorticity distribution. As long as the fluid velocity does not 
exceed the wave velocity, so there is no stagnation, the vorticity at a point only 
depends on the value of the stream function at that point. This dependency 
is described by what is known as the vorticity distribution. 


Lemma 2.2 (Vorticity distribution [4]). Suppose that u < c. Then there 
exists a function 7 such that iw = 7 ( 0 ) in 

A notable consequence of the existence of a vorticity distribution is that 
Equation (2.3a) is trivially satisfied, because 

{0, A0} = 0j;(-7(0))x - 0x(-7(V’))l/ = 0 

by the chain rule. Observe also that the condition in Lemma 2.2 is sufficient, 
but not necessary. By assuming the existence of a vorticity distribution, we 
will still obtain solutions of the water-wave problem, even if u < c is not 
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satisfied. In fact, the solutions that we will hnd can exhibit stagnation and 
critical layers. The introduction of the vorticity distribution is standard for 
rotational waves, and was used already in [9]. 

We shall consider the case where 7 is affine. By making a shift oi it is 
sufficient to consider the case of linear 7 . After scaling to unit depth and 
scaling away the gravitational acceleration, the stream function formulation 
(2.3) reduces to 

(2.4a) All) = ill 

(2.4b) + 7 = on 5, 

(2.4c) if) = mo on B, 

(2.4d) if) = mi on 5, 

where we have introduced the bottom B := {{x,y) y = 0} and the surface 
S := {{x,y) : y = 1 +r/(x)}. The parameter a in (2.4a) controls the vorticity, 
and will be assumed to be negative. For positive a, one-dimensional—but 
not higher-dimensional—bifurcation is possible. More discussion on this can 
be found in [ 11 ]. 

Observe now that the system (2.4) makes sense also in less regular function 
spaces than those specified in Proposition 2.1, and we will therefore allow for 
less regular (but still classical) solutions. More precisely, we will search for 
solutions 

7 G (M) and if) G 0 ^;^(H,), 

where /? G (0,1), and the subscript e signifies the subspace of functions which 
are even in the horizontal variable. The motivation for working in these 
Holder spaces is that Theorem 3.11 then holds. 

Remark 2.3 (Regularity). Due to Equation (2.4a) and elliptic regularity for 
the differential operator a — A, the stream function ■0 is analytic in 0^. In 
fact, we show in Theorem 2.5 that this is true even up to the boundary. 

2.3. Trivial solutions and flattening. The solutions of (2.4) that we shall 
construct will be small perturbations of steady flows that are parallel to the 
bottom. These parallel flows are the trivial solutions of (2.4), in the sense 
that 7 = 0 and the stream function 0 only depends on y. By integrating 
Equation (2.4a), we arrive at trivial solutions of the form 

(2.5) 0o(y, A) := 7 cos(|a |^/^(7 - 1)-h A), A = ( 7 , a, A) G M^, 

with corresponding Q(A), mo(A) and mi(A) determined from (2.4b)-(2.4d) 
as 

12 61 (j(Al - mo(A) = 7 cos(A - |a|^/2), 

2 ’ mi (A) = 7 cos(A). 

Our goal is to find nontrivial solutions of (2.4) for certain values of A, 
corresponding to these particular values of Q, mo and mi. Eor technical 
reasons which we will elucidate later in Remark 3.12, it is assumed that 

(2.7) 0Oi/(l) =-Ai|a|^'^^sin(A) 7 ^ 0. 

As in (2.7), we will often omit the dependence on A in our notation. 
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The main difficulty with the system (2.4) is that it is a free-boundary 
problem, which entails that the domain is a priori unknown. There are several 
ways of fixing the domain. Here, we will use the “naive” flattening transform 


G : (x,y) X, 




giving a bijection from the sets B and S onto 

14 = {(x, s) : s G [0,1]}, H := {(x, s) : s = 0} and S := {{x, s) : s = 1}, 

respectively. Using that rj G C'K|e(lK), we hnd that the map G is a G^’^- 
diffeomorphism, with inverse given by 

G“^(x, s) = (x, (1 r]{x))s). 

If we define ■0 on 14 by ■0 :=■)/; o G“^, then (2.4b) and (2.4a) become 


( 2 . 8 ) 


5 . - 

1 + ?? ^ 


0s 




= q; 0 in 14, 


+ 1] = Q on S, 


2(1 + 1 ?) 

in the new flattened variables, for which we have the following: 


Lemma 2.4 (Equivalence [11]). For functions rj G GK[e(lK) and if G GK,’e(14^), 
the stream function formulation (2.4) is equivalent to the transformed problem 
in (2.8) for rj G G^je (1^) and 

{0 e : 00=0 = mo, 00=1 = mi}. 

Moreover, in this setting, a pair {r],ip) = (O,0(s)) solves (2.8) if and only if 

0 = 00 - 


With the trivial solutions found and the flattening transform introduced, 
we now elaborate on Remark 2.3. Any solution which is sufficiently close 
to a trivial solution is in fact analytic, as long as (2.7) holds. The precise 
statement can be found in Theorem 2.5 below. 


Theorem 2.5 (Regularity). Suppose that a solution (i?, 0) of the problem 
(2.4) in G^(M) x G^(M) is such that the normal derivative dnif of the stream 
function vanishes at no point on the surface. Then we have the following: 

(i) The surface profile r] is analytic. 

(a) The stream function if extends to an analytic function on an open set 
containing Tl^. 

The assumption on dnif holds when A satisfies (2.7) and if is sufficiently 
close to 0o(-,A) in G^(14). 

Proof. We start by showing that g is analytic. For this, we will use the 
approach taken in [3], which is to apply [21, Theorem 3.2]. In [3] this was 
done under the assumption of no stagnation, but it is sufficient to assume 
that stagnation does not occur on the surface. This corresponds to the 
Shapiro-Lopatinskii condition for a certain elliptic system. 
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Let be the component of \ 5 that does not contain Proceed to 
dehne the function u: U 5 U —)• M by 

i ■= G S'uf];!-, 

- mi {x,y) e Qri'J S, 

and the differential operator L by L ;= a — A. Observe that Equations (2.4b) 
and (2.4d) imply that 

f{y, dni^) ■■= ^ {dni^f + y - 1 - Q = 0 

on S. All the assumptions of [21, Theorem 3.2] are now satished, with G '■= L 
and F{u) '■= Lu + ami (see the remark immediately after the theorem). We 
conclude that y is analytic. 

Note now that the differential operator L is strongly elliptic in the sense 
of [26, Equation (1.7)]. Equipped with the fact that y is analytic, we can use 
[26, Theorem A] to conclude that V’ extends to an analytic function on an 
open set containing 

The hnal part of the theorem follows because 

dn'ip = y/l+Jrffipyi-, y) = ^s(-, 1), 

where is bounded away from 0 as long as il) is sufficiently close to i/’o 

in C'^(ii), due to the assumption that Equation (2.7) holds. □ 

Remark 2.6. Theorem 2.5 is a local result at heart. It is clear from the proof 
that if dn^{xQ,y{xo)) ^ 0, then y is analytic in a neighborhood of xq. This, 
in turn, implies that extends analytically across the surface near the point 
(xo,r/(xo)). 

Remark 2.7. Recall that the stream function if: is analytic on 14^, regardless 
of whether the condition on dnip on the surface in Theorem 2.5 is satisfied. 
It is worth noting that this implies, through the implicit function theorem, 
that the streamlines are analytic curves away from stagnation points. 

2.4. The linearized problem. In order to linearize Equation (2.8) around 
a trivial solution i/iq, we write + (p, and introduce the spaces 

A = Ai X A 2 := (M) x G (A) : <^|,=o = ^\s=i = O} 

and 

F = X y2 := (M) X C'f,,(?^). 

We will write w = (y, (p) for elements of X. To capture our assumptions, it is 
convenient to dehne the sets 

O := {w G X : minr/ > — 1}, 

and, to enforce that a < 0 and (2.7) hold, 

47 := {(y, a. A) G : /U / 0, a < 0, 0 < A < tt} . 
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We now define the map T = : O x ^ y by 


(2.9a) J^i(ze,A) 


T2{w,A) 


{1 + r ] l ){ il>os + ( j > s ) 

2(l + r/)2 


+ rj — (5(A), 


dx ds 

l+r] 


(V'o + <^) + 


'4^0ss 

(1 + riY 


-a(V'o + <^), 


where Yo is as in Equation (2.5) and (5(A) is given in (2.6). In (2.9a), it is 
understood that the functions Yos and (pg are evaluated at s = 1. It is clear 
that is well defined and smooth as a map O x U ^ Y. We wish to solve 
the equation 


( 2 . 10 ) 


T{w, A) = 0. 


We obtain the linearized problem by taking the partial derivative of 
with respect to w at the point (0, A). This yields 

(2.11a) DwJ^i{0,A)w = Yos(f>s-Yos'n + 'n, 

DwJ^2i0, A)w = (A - a)0 - sYosVxx - ‘^'ipOssV, 

where it again is understood that the functions are evaluated at s = 1 in 
(2.11a). By introducing an isomorphism, in Proposition 2.8 below, we can 
transform into a simpler elliptic operator. For this purpose, define 

A2 := {</> G (P2;f (?^) ; = 0}, A := Xi X A 2 , 

where we have the inclusion X C A cY. We will typically use the letter (p 
for elements of A 2 . 


Proposition 2.8 (The T isomorphism [11]). The bounded linear operator 
T (A) : A 2 —)• A defined by 




fi\s=^ 


is an isomorphism of Banach spaces, and the operator 


£(A) = (A(A), £ 2 (A)) := Zl^A(0,A)r(A) : A 2 ^ A 


satisfies 

( 2 . 12 ) 


C{X)cP 




, (A - a)(p 


S=1 


Proof. That T is well defined and an isomorphism is almost immediate. The 
expression for £(A) in (2.12) follows by direct computation. □ 


3. The kernel and dimensional reduction 


Introduce the complex parameter 


e 


n 


6{n, a) := \/ a + A A = 


i^/\a\~--nf~^Y 


n > 
n < 
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for nonnegative integers n. This parameter will appear in functions of the 
form cosh(0„s) and sinh(0„s)/0„, which are always real-valued. We record 
that 


cosh( 0 „s) 

sinh( 0 „s) 


COsh( 0n s), 

n > 

a 

COs(|0n|s), 

n < 

a 

sinh(|6'n|s)/|6'n|, 

n > 

a 

sin(|6'n|s)/|6'n|, 

n < 

a K. 


In the event that 6n = 0, we will interpret expressions with 6n as extended 
by continuity. In particular, sinh(0„s)/0n is interpreted as s. 

We now describe the kernel of >C(A), which is directly related to the kernel 
of A) through T(A), in terms of the above functions. The following 

proposition is stated, but not proved, in [11]. We include its proof because it 
is instructive. 


Proposition 3.1 (Kernel of /2(A) [11]). Let K ^U. A basis for ker/2(A) is 
then given by {4>n}n£M, where 


(3.1) 


4>n{x, s) := cos(nKx) 


sinh(0,i 


and M = M (A) is the finite set of all n G No satisfying the kernel equation 
(3.2) l{n,a)=r{A), 


where 


(3.3) 

Proof. Suppose 


fin, a) ■■= 6nCot\i{6n), 
r(A) := ^ 


-|- cot(A). 


tif\oi\ sin^(A) 

that (j) G ker/2(A), and expand it in a Fourier series 


OO 

(t){x,s) = E an{s) cofinKx). 

n=0 

From C 2 {A)(p = 0, we deduce that the coefficients satisfy 
(3.4) af{s)-0nan{s) = O, sG(0, 1), 

while 4>\s=o = 0 and Ci{A)(f = 0 yield the boundary conditions 
(3.5a) an(0) = 0, 

(3.5b) - ^V’Oss(l) 

for all n > 0 . 

The general solution of (3.4) with the boundary condition (3.5a) is 

, , „ sinh(0„s) „ _ 

an{s) = Bn - n -> Bn G M, n > 0, 

for which the Robin condition (3.5b) reduces to 

cosh(»„) - (*„„(!) + 55^) B„ = 0. 
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Hence, if Bn (and thus a^) is nonzero, then 

(3.6) V’Os(l) cosh(6»„) - ^V’Oss(l) + sinh(6*n) ^ ^ 

must hold. Observe that Equation (3.6) implies that smh{9n)/0n 7^ 0; 
otherwise we would have cosh(0„) = sinh(0„) = 0, and therefore exp(0jj) = 0. 
Thus, by inserting the definition (2.5) of into Equation (3.6), we arrive at 
(3.2). This condition is also sufficient for (pn to lie in the kernel. 

The set M of n G Nq such that (3.2) holds is hnite, because the function 
l{-,a) is strictly increasing as soon as n > /k. □ 

Remark 3.2. In order to get nontrivial solutions, A should be chosen such 
that 0 ^ M{K). The function (pQ, see (3.1), does not depend on x. 

The next lemma, inspired by [19, Theorem IV.5.17], serves to show that 
the set-valued map M: U ^ 2^° defined in Proposition 3.1 is upper semicon- 
tinuous. This implies that no new solutions of the kernel equation (3.2) can 
appear if A is perturbed slightly. 

Lemma 3.3 (Upper semicontinuity). Let A* G U. Then 

M(A) c M{A*) 

for all A in a neighborhood of A*. 

Proof Suppose that this is not the case. Then we can construct a sequence 
(Aj)jgf^ converging to A*, and a corresponding sequence (ni)igN such that 
Ui ^ M{A*) and /(n*, ai) = r(Aj) for all i G N. By the continuity of r at A*, 
the sequence (r(Aj))jgf^, and therefore {l{ni,ai))i£}^, is bounded. This implies 
that (ni)jgN is bounded, so we may assume that it is constant. Thus there is 
an n ^ M{A*) such that l(n, ai) = r{Ai) for all i G N. The boundedness of 
the sequences now ensures that l{n,-) is well-defined and continuous at a*. 
We conclude that l{n, a*) = r(A*), which contradicts n ^ M{A*). □ 

T 

h3{a,X) 
h2{a,X) 

Hi{a,X) 


a 

A 



Figure 1. An illustration of Theorem 3.4 when |M(A*)| = 3. 

A full description of the kernel can be given near A*. 

We can now use Lemma 3.3 to give a local description of the structure of 
the kernel equation. This will be useful when describing the solution set of 
(2.10). See also Figure 1. 
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Theorem 3.4 (Local description). Suppose that M(A*) 
Then we may define 

(3.7) ft(a,A):=—^- 

sin(A)(/(ni, a) — cot(A))^/2 

on a neighborhood of (a*, A*), 

^*(A) ;=A) (=••• = A)) 

on a neighborhood of X*, and we have 


M(A) 


' M{A*) 
0 

\ 


for all K in a neighborhood of A*. 


a = a*, pL = pifiX), 
a a*, p, = p,i{a, X), 

otherwise, 


{ni,... ,nAr}. 
I <i<N, 


Proof. By Lemma 3.3 there is a neighborhood of A* in which M(A) is the 
set of Ui £ M{A*) for which l{ni, a) = r(A). Observe now that A sufficiently 
close to A* we have /(n*, a) = r(A) if and only ii p = pfia, A), where pi is as 
in (3.7). Moreover, if i j then l{ni, ■) — l{nj, ■) is a nonzero analytic function 
on a neighborhood of a*. It follows that we may choose the neighborhood 
of A* in such a way that the only intersection of the graphs of the pi occurs 
when a = a*. □ 


The bifurcation results in Sections 4 and 5 are valid under the assumption 
that ker£(A) is respectively one- and two-dimensional. Lemma 3.5 below is 
a general result on the kernel equation (3.2) from [11], which in particular 
shows that it is indeed possible to choose A £hl such that the dimension of 
the kernel is one or two. 


Lemma 3.5 (Kernel equation [11]). 

(i) For every a and any n for which fin, a) is well-defined there are p and 
X such that n G M(A). 

(a) Suppose that X G [7r/2,7r)an4 that ni,n 2 G No satisfy 

(£) 


Then there are a and p such that ni,n 2 G M(A) and any other solution 
of (3.2) must be smaller than ni. 


It is, however, the case that higher-dimensional kernels are, in a sense, 
rare^: 


Lemma 3.6. Let J be set of all values of a for which there exist X and p 
such that |M(A)| > 2. Then the limit points of J are contained in the set 

(3.8) { —(mf -|- m^vr^) : mi G No, m 2 G N}. 

In particular, J consists of isolated points, except possibly those that lie in 
the set defined in (3.8), and has countable closure. 

variation of this was pointed out already in [11]; however, not taking into account 
the points where I is not well-defined. We slightly improve upon the result here. Both J 
and its closure are small. 
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Proof. The set defined in (3.8) consists precisely of the values of a for which 
there is at least one n G No such that l(n,a) is not well-defined. Let a be 
such that it is not in this set. 

Suppose, to the contrary, that there is a sequence (ai)igis} satisfying a* 7 ^ a 
for all i G N and which converges to a, with corresponding sequences 
and (n 2 ,i)ieN such that ni^i < n 2 ,i and 

for all i G N. We must necessarily have 

Rl,i < 

for all i G N, and so the sequence is bounded. The continuity of 

/(n, •) at a for each n G Nq now implies that (/(ni^i, Q;j))igN, and therefore also 
(^(n 2 ,i, cii))ieN) is bounded. This, in turn, implies that (n 2 ,i)ieN is bounded. 

By going to a subsequence, we may assume that both (ni,i)jgN and (n 2 ,i)ieN 
are constant. Thus there are m < n 2 G No such that l{ni, a) = l{n 2 , a) and 

l{ni,ai) = l{n2,ai) 

for all i G N. But this is impossible, because l{n 2 , ■) •) is a nonconstant 

analytic function in a neighborhood of a. □ 

For later use, we give some explicit examples of one- and two-dimensional 
kernels for /2(A). All satisfy r(A) = 1, and the two-dimensional examples 
have been chosen such that 9{n2, ot) = 0. To simplify the parameters involved, 
we choose specific values of k. 

Example 3.7 (Explicit kernels). Let a be the smallest positive solution of 
the equation xcot(x) = 1 . 

(i) When k = 1, ^ = 1, a = —1 and A = 7 r/ 2 , the kernel is one-dimensional, 

being spanned by s) = cos(x)s. 

(ii) Let K = When /j, = 1/{2 k), a = —4p and A = 7 r/ 2 , the kernel 

is two-dimensional, with M = {1,2}. 

(iii) Let k = When ji = l/(3/t), a = —and A = 7 r/ 2 , the kernel 

is two-dimensional, with M = {2,3}. 

3.1. Arbitrarily large kernels. We now address a question that was raised 
in [11]: Do there exist A G ^ such that kerA) is at least three- 
dimensional? By also letting the wave number k vary, this question was 
answered in the affirmative for dimension three in [12]. In essence, their 
result says that many two-dimensional kernels can be modified in order to 
yield a three-dimensional kernel. Here, we use a different approach to find 
kernels of arbitrary dimension for any k in a set K that is dense in (0, 00 ). 

For any a < 0 and A G p/2, tt), we can obtain r(A) = 0 by choosing /x to 
satisfy 

2 _ 2 

sin(2A) 

which reduces the kernel equation (3.2) to finding m G N and n G Nq such 
that 
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is satisfied. This can be written in the form 

(3.9) (2n-y+ {2m-lf = ^. 

V TT/ TT^ 

We first consider the case k = n. 

Lemma 3.8 (Arbitrary kernel with k = ir). For k = n and any N £ N, 
there exist K £14 sueh that |M(A)| = N and 0 ^ M(A). 

Proof. When k = tt and a = for an odd number H £ N, (3.9) 

becomes the Diophantine equation 

(3.10) {2nf + i2m-lf = H. 

The size of the kernel then corresponds to the number of representations of 
H as the sum of two squares. As long as H is not a square number, any such 
representation has n ^ 0 (see Remark 3.2). 

In order to conclude, we therefore need to find an odd non-square number 
H such that H has exactly N representations as a sum of squares. By [15, 
Theorem 3 in Chapter 2], this is the case for instance when H = fQj. 

prime p £ 4N -|- 1. □ 

Some examples of Lemma 3.8 are listed below, for various choices of H 
in (3.10). The values of H used in second and third example, which are the 
smallest possible, are not in the form They can easily be deduced by 

using the general formula given in [15]. 

Example 3.9 (Higher-dimensional kernels). 

(i) The choice H = 5^'^“^ = 3125 yields a three-dimensional kernel, with 
M = {5,19, 25}. However, this is not the smallest example: 

(ii) Since 325 = 6^ -|- 17^ = 10^ -I- 15^ = 18^ -|- 1^ (with no other represen¬ 
tations), the choice H = 325 yields a three-dimensional kernel, with 
M = {3,5,9}. 

(hi) Since 1105 = 4^ + 33^ = 12^ 31^ = 24^ -h 23^ = 32^ 9^, the choice 

H = 1105 yields a four-dimensional kernel, with M = {2,6,12,16}. 

Let denote the set of positive rational numbers with odd numerators 
when reduced to lowest terms. We can then generalize Lemma 3.8 in the 
following way: 

Theorem 3.10 (Arbitrary kernel). For k £ ttQq and any N £'N, there exist 
K £lA sueh that |M(A)| = N and 0 ^ M(A). 

Proof Write k = nrfs, with r and s coprime. When a = —7r^r^77/4, (3.9) 
becomes 

(3.11) r^(2n)^-|-s^(2m — 1)^ = r^s^i7. 

Choosing H = £qj. pj-ij^e p £ 4N -|- 1, we know that (3.10) has exactly 

N solutions {mj,hj) in N^. The pairs {mj,nj) £ defined by 

2mj — 1 = r{2rhj — 1), nj = shj 

then solve Equation (3.11). 

Moreover, these are the only solutions: Suppose that (m, n) solves Equa¬ 
tion (3.11). Then r | {2m — 1) and s \ 2n, by coprimality of r and s. 
It follows that 2m — 1 = r{2m — 1) and 2n = sh, where m and h solve 
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n? + (2m — 1)^ = H. Since H is odd, n = 2h. Uniqueness in Equation (3.10) 
now yields the result. □ 

Although we provide kernels of arbitrary dimension in Theorem 3.10, the 
corresponding triples A satisfy r(A) = 0, unlike the three-dimensional kernels 
obtained in [12]. In particular, this means that the two-dimensional bifurca¬ 
tion result in Theorem 5.2 does not apply for the kernels from Theorem 3.10 
with N = 2. We remark that an obstacle for higher-dimensional bifurcation 
is that there are only four parameters to work with, namely A and k. This 
may be remedied by for instance including surface tension. 

An application of Theorem 3.10 is one-dimensional bifurcation with several 
different wave numbers for fixed A. If the set M(A) = {ni,..., njsf} is such 
rii f rij for all i ^ j, we can make restrictions to each and then apply 

Theorem 4.3. This will yield N different solution curves. Two examples 
for which the condition on M is fulfilled are H = 725 and H = 3145, 
corresponding to M = {5, 7,13} and M = {18, 24, 26, 28}, respectively. 

3.2. Lyapunov—Schmidt reduction. Before we can reduce (2.10) to a 
finite-dimensional problem by applying the Lyapunov-Schmidt reduction, we 
need the following result from elliptic theory. 

Theorem 3.11 (Fredholm property [11]). The operator C{A) is Fredholm for 
each A £U, with index 0. The range of C{A) is the orthogonal complement of 

Z '■= { (f) : (f ^ ker £(A)} C X C T 

in Y with respect to the inner product 

‘I'K j K 1 I'K j K 

(3.12) {wi,W2)y = Jriir]2dx + J J 4)i(f>2dxds, Wj^Y. 

0 0 0 

Let Wn ■= (r/ 0 ^,(/>„) for n G M(A), where 4>n and M(A) are as in Proposi¬ 
tion 3.1. Then the projection Hz :Y^Z onto Z along ran£(A) is given 
by 

fo,o\ rr 'iW {w,Wn)Y ~ 

(3.13) Hzw = 2^ , 2 

n£M{A) 

Remark 3.12. That I1^J^(0, A) be Fredholm is the main reason for making 
the assumption (2.7). When we have equality in (2.7), (2.11a) reduces to 
Du,jFi{0, A)w = rj, whence ranL)^t,J^i( 0 , A) = Xi. The operator Z1^J^(0, A) 
then cannot be Fredholm, since Xi is not closed in Yi. 

Let A* G ^ be a triple {fi* , a*, X*) such that N := |M(A*)| > 1. Then 
Proposition 3.1 says that the pairs 

< := r(A*)</.; G X, neM{A*), 

span the kernel of DwiF{0, A*). Since the kernel is hnite-dimensional, there 
exists a closed subspace Xq C X such that 

X = ker71^X(0,A*)©Xo. 

By Theorem 3.11, we can also decompose Y into the direct sum 

Y = Z Oran C{A*), 
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where Z := span{A*), since these are orthogonal complements in 
the inner product (3.12) on Y. Applying the Lyapunov-Schmidt reduction 
(see e.g. Kielhofer [20]) for these decompositions of X and Y, we obtain the 
following lemma. 

Lemma 3.13 (Lyapunov-Schmidt). There exist open neighborhoods N of 
0 in kei DwT'{0, A*), Ai of 0 in Xq, and li' of A* in U, and a uniquely 
determined function if : N x lA' ^ Ai such that 

A) = 0 for w £ M + Ai, A £ U', 

if and only if w = w* + 'f{w* ,A) and w* = YlneM{A*) ^nWn ^ solves the 
finite-dimensional problem 

d>(t,A) = 0 fort £V, A£U', 

where 

^{t,A):=IlzT'{w,A) and V := {(tn)neM(A*) £ : w'* £ AA}. 

The function is smooth, and satisfies '0(0, A) = 0 for all A £ W, and 
Z1^0(O,A*) = 0. 


4. OnE-DIMENSIONAL BIFURCATION 

We are now in a position to show that a curve of nontrivial solutions 
of (2.10) bifurcates from each point (0, A*) £ X xU where the kernel of 
Du,T'{0, A*) is one-dimensional, given that A* satisfies an additional technical 
condition. This condition comes from Lemma 4.1 below. 


Lemma 4.1 (Orthogonality). Suppose that n £ M{A), so that the function 0„ 
given by (3.1) lies in ker£(A). Then, if Wn '■= is the corresponding 

basis function of Z, we have 


(4.1) 

where 


{DxT[X)(f)Yi,'Wn)Y — ^ 


sinh(6'n) 


A := -- 


2tt 


does not depend on n. In particular. 


cot(A) -|- 


2m,|3/2 


pL \a 


{DxC{A)(pri,Wn)Y = 0 */ COt(A) = 


2 


Proof. Recalling (2.12), we find the derivative 

Z1a/:(A)0= (^0o.a(1)0.|.=i- (0Os.a(1)-||^)0|s=i,o). 
Using that (pn{x, s) = cos(nKx) sinh(0„s)/0n, we get 

DxCi{A)4)n = cos(nKx), 
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where 


A ■■= 'tpos\(i)l(n,a) - ^V'ossa(I) - ^ 


cot (A) + 


/l2|a|3/2 


V’Os(l) 

by the kernel equation (3.2) and the definition of t/^o- 
Since Wn = (r/^„,(pn), with 


Wn(x) = - 


we now find 


4>n{x, 1) _ 1 sinh(6<n) 

'IpOsi'i-) V’Os(l) 9n 

21^ j K 


cos(nKx), 


{DxC{A)(pn,Wn)Y = j 'n<f>„D\Ci{A)4>n dx 


ttA / sinh(6*n) 

K1p0s{l) V 


which is (4.1) with A = —TrA/{K'il>osi'^))- 


□ 


We will refer to 

(4.2) cot(A)74-'ddfil 

as the transversality condition, because it corresponds to transversality in the 
Crandall-Rabinowitz theorem (see [8] or [20]). Note that all the examples we 
provided in Example 3.7 satisfy this condition. It is straightforward to check 
that the transversality condition fails at A* gU precisely when /r(,(A*) = 0 
in Theorem 3.4. This means that we can obtain the following by moving 
slightly along the graph of /U* (see Figure 1). 


Lemma 4.2. Suppose that A* = {p,*,a*, X*) gU is such that the transver¬ 
sality condition (4.2) fails. Then there are p, X gM. with A = {p, a*, X) G U 
such that the transversality condition holds and M(A) = M(A*). The triple 
A can he chosen arbitrarily close to A*. 


The one-dimensional bifurcation result is an application of the Crandall- 
Rabinowitz bifurcation theorem. To clarify the proof of the two-dimensional 
bifurcation in the next section, we will nonetheless spell out the details of 
the proof. 

Theorem 4.3 (One-dimensional bifurcation). Suppose that A* gU is such 
that M{A*) = {n} with n G N, and therefore that 

kei DwJ^{0, A*) = span {re*}, 

where w* = T{A*)(p*, with (p* := cpn o,s in Proposition 3.1. If the transversality 
condition (4.2) holds, there exists a smooth curve {(w{t), X{t)) ; 0 < |t| < e} 
of nontrivial small-amplitude solutions to 

(4.3) T{w,p*,a*,X) = 0, 
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in O X (0, tt), passing through (u;(0),A(0)) = (0,A*), with 

(4.4) w{t) = tw* + 0{t^) in X as t ^ 0. 

These are all the nontrivial solutions of (4.3) in a neighborhood of (0, A*) in 
O X (0, tt). 

Proof. Using Lemma 3.13, we know that there exists a neighborhood of 
(0, A*) in O X (0, tt) for which the equation X{w, //*, a*, A) = 0 is equivalent 
to 4>(t, n*,a*,X) = 0, where t G M. From the same lemma we also have the 
identity <h(0, A) = 0, and hence we can write 

4>(t,A)= [ dz{^{tz,A))dz = t'^{t,A), 

Jo 

where 

(4.5) 'I'(t, A) := f ^t{tz,A)dz 

Jo 

is smooth. For nontrivial solutions {t ^ 0), the equations = 0 and 4' = 0 are 
equivalent, whence we need only concern ourselves with the latter equation. 

We want to apply the implicit function theorem to ik, which requires that 
4^(0, A*) = 0 and 4';^(0, A*) ^ 0 (recall that Z is one-dimensional). Now, 
from (4.5), we find 

4'(0,A*) = $j(0,A*), 

4/a( 0,A*) = 4>a(0,A*), 

so these are the derivatives of we need to compute. By the definition of <k, 

(4.6) d>t{t, A) = UzD^X{tw* + i;{tw*,A),A){w* + D^f;{tw*,A)w*), 

and so by evaluating in t = 0, and using the properties of fr listed in 
Lemma 3.13, we have 

^t(0, A) = K){w* + D^TfiO, A)w*), 

which also yields 

4>u(0, A) = UzD^xJ^iO, X){w* + D^fT{tw\K)) 

+ \dzDy,F{d, A)Dy,x'4){w*,A). 

We now obtain 

^(0, A*) = $t(0. A*) = IlzD^X{d,A*)w* = 0, 

because Dw'il>(0, A*) = 0 by the last part of Lemma 3.13, and because 11^ 
projects along T:anDwX{0, A*). Similarly, 

(4.7) ^a(0, A*) = UzD^xJ^iO, A*)w*. 

Note that Diujz{0, A)w* = £,{A)T{A)~^w*, and hence 

D^xJ^{0,A*)w* = DxC{A*)cj)* - D^X{d,A*)dxT{A*)ctT\ 
which implies that Equation (4.7) can be written 

^tx{d,A*) = YizDxC{A*)(f*, 



18 


AILO AASEN AND KRISTOFFER VARHOLM 


again using that Tlz projects along ran A*). We can now use 

Lemma 4.1 to deduce that 'La(0, A*) = <^^^(0, A*) ^ 0, due to the assumption 
of transversality. 

Finally, since 'I'(0,A*) = 0 and 'La(0,A*) ^ 0, we can invoke the im¬ 
plicit function theorem to deduce that there exists an e > 0 and a smooth 
function A; {—s,e) —)• M with A(0) = A* such that fj,*,a*, X{t)) = 
0. Moreover, the curve {{t,X(t)) : |t| < e} describes all solutions to 
fi*,a*, X) = 0 in a neighborhood of (0,A*). The corresponding solu¬ 
tion curve to F{w, fi*,a*,X) = 0 is {(w{t),X{t)) : |t| < e}, where w{t) := 
tw* -)- , fi*, a*, X{t)). It follows that 

w(t) = w* + Dwipitw*, fi*, a*, X{t))w* + , fi*, a*, X{t))X{t), 

and we can conclude, once again using the properties of given in Lemma 3.13, 
that w{0) = 0 and w{0) = w*. Consequently, we obtain (4.4). □ 

If e is sufficiently small, the waves obtained from Theorem 4.3 are Stokes 
waves. This can be seen from the asymptotic formula in (4.4). 

4.1. Properties of the bifurcation curve. The one-dimensional bifurca¬ 
tion result in Theorem 4.3 is analogous to [11, Theorem 4.6], which uses fj, 
instead of A as the bifurcation parameter. Other than the parameters, the 
main difference between the theorems is the addition of the transversality 
condition (4.2) for bifurcation with respect to A. Here, we will investigate 
the properties of the solution curves more closely. 

The motivation is to understand the solution set of (2.10) better, and 
in particular to rule out the possibility that the solution curve found here 
coincides with the one from [11]. The only way this can occur is if X{t) and 
in the notation of [11], are constant along the curves. (If they were 
constant, we would obtain the same solutions by uniqueness in Theorem 4.3.) 
Proposition 4.4 shows that we need to consider at least second-order properties 
of the bifurcation curve in order to achieve this. 

Proposition 4.4 (First derivative of A). Under the hypothesis of Theorem 4-3, 
the function X satisfies 

A(0) = 0, 

and so the bifurcation parameter is constant to the first order along the 
bifurcation curve. 

Proof We adopt the notation used in the proof of Theorem 4.3. Differentia¬ 
tion of the identity 'I'(t, p,*,a*, X{t)) = 0, and evaluation at t = 0, yields the 
equation 

(4.8) Ti(0,A*) + ^A(0,A*)A(0) = 0 

for the derivative of A at the origin. From (4.7) and the discussion immediately 
after, we know that 

(4.9) Ta(0, a*) = UzDxC{A*)f* / 0, 

which means that (4.8) uniquely determines A(0). However, we still need to 
compute Tt(0, A*). 
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From (4.5), we obtain 

and differentiation in Equation (4.6) leads to 

A) = UzDlT{tw* + ^P{tw*,A),A){w* + A)w*f 

+ YlzD^Fitw* + ^{tw*,A), A)Dli){tw\A){w*f. 


Hence, by using the properties of V’ given in Lemma 3.13, and using that liz 
projects along the range of DyjT{Q, A*)^ we hnd 

(4.10) 4-^(0, A*) = hizDlT{Q, A*){w*f. 

Using Equations (4.8) to (4.10) and the formula (3.13) for Hz given in 
Theorem 3.11, we now hnd 


(4.11) 


{DlmA*){w*)\w*)y 

2{DxC{A*)4>*,w*)y 


and so it is sufficient to show that the numerator. 


(4.12) 


27r//t 1 27 t/k 

[ r]*DlTi{0,A){w*fdx+ [ [ 4>*DlT2{0,A){w*fdxds, 


0 0 


vanishes. Since w* = T{A*)(j)* with cp* being a separable function of 
X and s, so is t]* . Moreover, we see from Equation (3.1) that their x- 
dependence is through cos(n«:x). Thus each term in A*)(w*)‘^ and 

A*){w*)‘^ has an x-dependence of the form sin“(nAx) cos^(nKx) with 
a + 6 = 2 (see the derivatives listed in Appendix A). It follows that we will be 
integrating terms whose x-dependence is sin“(nKx) cos^(nKx) with a + 6 = 3 
in (4.12), and therefore that the numerator in (4.11) vanishes. □ 


Remark 4.5. Equation (4.11) also holds if /x is substituted for A. Since the 
proof of Proposition 4.4 only depends on the fact that the numerator in (4.11) 
vanishes, we can conclude that one also has /l(0) = 0 when using fi as the 
bifurcation parameter. 


To consider the question of second-order behavior of the solution curve, let 
us return to the expression for w we found in the one-dimensional bifurcation 
result Theorem 4.3, namely 

(4.13) w{t) = w* + Dw^w* + 

where Duj'ip and tpx are evaluated at {tw*, n*, a*, X{t)). Taking another 
derivative in (4.13) yields 

(4.14) w{t) = + 2X{t)Dwx'ipw* + X{tfipxx + X{t)^px- 

This simplifies signihcantly at t = 0. An expression for A(0) can also be found, 
akin to how (4.11) was derived. The details are omitted, see for instance [20, 
Section 1.6]. 
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Proposition 4.6 (Second derivatives). Under the hypothesis of Theorem 4-3, 
we have 

W{0) = 

and 

3{DxC{A*)(j)*,w*)Y 

for the solution curve {w{t), X{t)). 

Proof for w{0). In view of Lemma 3.13, we have that '0(0, A) = 0 for all A in 
an open neighborhood of A*, and thus 0;^(O, A*) is zero. Finally, the second 
and third terms in (4.14) vanish due to Proposition 4.4. □ 


Remark 4.7. The expression for can be obtained by simply substituting 
fi for A in the expression for A(0). From this, it follows that 


A(0) 


{D\C{A*)(j)* ,w*)y ^(cot(A) + /2) 


In particular, this implies that A(0) and /u(0) must have either the same or 
opposite sign, depending on the sign of jj, and which “side” of the transversality 
condition (4.2) A* is on. 


We now give a more transparent description of 11^0(0, A*)(rc*)^, which 
Proposition 4.6 shows is required for computing both u;(0) and A(0). 


Lemma 4.8 (Description of iA^0(O, A*)(t(;*)^). Write 

iA^J'i(0, A*){w*)‘^ = Co + C 2 cos{2nKx), 

(4.15) „ „ 

01 ^ 7 ^ 2 ( 0 , A*){w*)^ = 60 (s) + 62 (s) cos(2nKx), 


‘2, 8 ^ 

and let f G C^’e (D) he such that 

C{x, s) := ao{s) + 02 ( 5 ) cos(2nKx), 
where the coefficients oq and 02 solve the boundary value problems 

a'jis) - 0]naj{s) = -bjis), 

0Os(l)a' (l) - ^0oss(l) + ^^^^^^ ^ai(l) = -Cj, 

Dl^f{0,A*){w*f = T{A*)(:. 

Proof. That D‘^iF{0, A*)(w*)‘^ can always be written as in (4.15) can be 
deduced from the expressions for the derivatives of T listed in Appendix A. 
The function 0 satishes the identity 

(I - Iiz)F{tw* + f;{tw*,A*),A*) = 0 

for sufficiently small t. If we take two derivatives of this equation and evaluate 
at t = 0 we obtain the equation 

(/ - Uz){DlT{0, A*){w*f + Dy,T{0, A*)Dli;{0, A*)iw*f) = 0 


(4.16) 

a,(0) = 0, 
for j = 0,2. Then 
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for Since we established in the proof of Proposition 4.4 

that lies in the range of A*), and since 11^ projects 

along ranA*), this implies that 

(4.17) A*)Dli;{0, A*)iw*f = -DIt{0, A*)(w*f, 

which uniquely determines £>^1^(0, A*)(w*)^. 

If we now introduce the function ^ G X 2 by 

T(A*K:= 

then Equation (4.17) can be written 

£(A*)C = -DlF{Q,A*){w*f. 

Utilizing (2.12), this proves the lemma. □ 

Due to the form of w* and the expressions for the derivatives of J- in 
Appendix A, we know that the coefficients uq and 02 in Lemma 4.8 are 
polynomials in s, smh(9jns)/9jn and cosh(0jjjs) for j = 0,1,2. They can, 
with some effort, be computed explicitly using a computer algebra system. 
However, the general expressions are much too long to perform any useful 
analysis of the second derivatives. We will therefore content ourselves with 
presenting the result for the hrst special case of Example 3.7, which was 
constructed specihcahy to make (j)* and 'ipQ as simple as possible. This, in 
turn, yields particularly simple oq and 02 . 


Theorem 4.9 (Special case). When k = 1 and A* = (1,-1,7r/2), the 
functions oq and 02 are given by 


(4.18) 


, , 1 o . , _ 3sin(s) 

ao{s) = s + -s sm(s - 1 ) + - --r- ^ 7 -^, 

2 2 (cos(l) — sm(l)) 

In./ ..x sinh(\/3s) 

= » + -» - 1) + 2(V3cosh(G!)-sinh(V!))’ 


respectively. This yields 

(4.19) A(0) — 2 3ao(l) + 2 ^ 2 ( 1 ) < 0- 


In particular, A is not constant along the bifurcation curve, which therefore 
does not coincide with the one found in [11]. 


Proof. The kernel of £(A*) is spanned by cos(x)s, and moreover V’o('S) = 
— sin(s — 1). One may check that for this special case, the coefficients in 
Equation (4.15) are given by 

Co = 2, bo{s) = —s — 2scos(s — 1) — sin(s — 1), 

C 2 = 1, ^ 2 ( 3 ) = 3s — 2s cos(s — 1) + (2s^ — 1) sin(s — 1). 

It follows by direct verihcation that the functions oq and 02 in (4.18) solve the 
boundary value problems in (4.16). Finally, a long (but direct) computation 

from the expression for A(0) in Proposition 4.6 yields (4.19). □ 


For the same special case as in Theorem 4.9, we can consider bifurcation 
from other points on the graph of the associated function that was 
introduced in Theorem 3.4. One may verify that the numerator in the 
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Figure 2. The solution curves emanating from the graph of 
making up the solution set of T = 0 when a is fixed. The 
specihc point used in Theorem 4.9 can be found to the left of 
where the transversality condition fails. 

expression for A(0) in Proposition 4.6 is negative on the entire graph of ;U*. 
This means that, locally, the solution set of the function T from the proof of 
Theorem 4.3 looks qualitatively like the surface shown in Figure 2, when a is 
fixed. Recall that the transversality condition corresponds to jj,'^ ^ 0, and 

observe that A(0) changes sign when this condition fails. 

We remark that for some other choices of k and A the numerator does 
change sign on the graph of It follows that Figure 2 does not, in general, 
tell the whole story. 

4.2. Local description. Using [11, Theorem 4.6], or Theorem 4.3 when 
the transversality condition (4.2) is fulfilled, we can describe all solutions 
of Equation (2.10) in a neighborhood of any (0,A*) in X xli for which 
\M{A*)\ = 1. 

Suppose that we have such a point, and that M(A*) = {n}. Then Theo¬ 
rem 3.4 tells us that there is a neighborhood of A* in which 

M(A) = |'^”^ M = Mi(«,A), 

10 otherwise. 

This allows us to invoke [11, Theorem 4.6] on each point on the graph of //i, 
obtaining a family of solution curves. These are, in fact, all the nontrivial 
solutions near (0, A*): 

Theorem 4.10 (Local description). The above family S of solution curves bi¬ 
furcating from points (0, ^i(a, A), a. A) for (a. A) in a neighborhood of {a*, A*) 
contains all nontrivial solutions of Equation (2.10) in a neighborhood of 
(0,A*) inXxU. 

Proof For each A = (//i(a. A), a. A) we have uniqueness in a set 

U (A) ;= {(re, p!,a,X) ^ X xU : ||rc|| V — /Ui(a, A)| < 5(A)}, 

in the sense that all nontrivial solutions in U{X) are given by the solution 
curve obtained in [11, Theorem 4.6]. Due to the regularity of the problem, 
and by possibly shrinking the neighborhood of (a*, A*), we can assume that 
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/il(a, A) 


Figure 3. All nontrivial solntions near (0, A*) can be fonnd 
by bifurcation from points on the graph of 


(5(A) is constant. It is then clear that the family S yields all the nontrivial 
solutions in the open neighborhood 

[/ := U ^(A) 

of (0,A*), see Figure 3. □ 

Remark 4.11. We mentioned above that the same procedure can be performed 
using Theorem 4.3 instead of [11, Theorem 4.6] when the transversality 
condition is fnlfilled. The implication is that, locally, the same solntions can 
be fonnd throngh bifnrcation with either // or A. It is not clear whether this 
is still the case for possible global solution curves. 


5. Two-dimensional bifurcation 


For two-dimensional bifurcation we will use a as the second bifurcation 
parameter. We therefore need the following analogue of Lemma 4.1 for the 
parameter a. 


Lemma 5.1. Suppose that (pn £ ker£(A), where (pn is as defined in (3.1). 
Then, if Wn = iv<l>nj4’n) is the corresponding basis function of Z, we have 

(Zl„£(A),/.„,u)„)y = +/(0„,), 


where 


and f is defined by 



1 

sin^(A) 


cot(A) 

21^. 


fit) ■■= 


TT 

K 


ft — cosh(t) sinh(t) 



2t3 


t 0, 
t = 0. 


Suppose that M{A*) = {ni,n 2 }, where ni < n 2 , and that A* satisfies 
the transversality condition (4.2) (which also appears for two-dimensional 
bifurcation). Denote the subspace of X consisting of functions that have 
wavenumber nn in the horizontal variable by Then ker £(A*)|= 

span{(/>* 2 }- Using the local description from Section 4.2, we obtain the 
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set of all nontrivial solutions of (2.10) in a neighborhood of (0,A*) 

in xM X L/. Similarly, under the condition that ni f n 2 , the kernel of 
£(A)|j^(nj) is spanned by (/>^^, and we obtain the set of all nontrivial 

solutions in a neighborhood of (0, A*) in xU. The next result describes 

bimodal solutions near (0, A*), which are neither in nor in 

As mentioned before Lemma 5.1, we will use a as the second bifurcation 
parameter. Hence, we will look for solutions of the equation 

(5.1) /r*, a. A) = 0 

ioT w £ O and {ji*,a,\) G U. Let therefore, for j G {1,2}, the set 
consist of all {w,a,\) such that {w, fi*,a, X) G 

Theorem 5.2 (Two-dimensional bifurcation). Suppose that A* £U is such 
that the transversality condition (4.2) holds, and that 

kerll^j^(0. A*) = span{rt;{,rc2}, w* = T{A*)(j)*, 

with 1 < ni < n 2 and 0} := (j)^^ as in Proposition 3.1. Furthermore, suppose 
that either r(A*) ^ {0,1} or 0{n2,a*) = 0 (in which case r(A*) = 1). 

(i) If ni \ n 2 , there exists a smooth family of nontrivial small-amplitude 
solutions 

Sf,* := {{w{ti,t2),a{ti,t2),X{ti,t2)) : 0 < |(ti,t2)| < e} 

of ( 5 . 1 ) in O X (—00,0) X (0, tt), passing through (0,a*,A*) when 
(H, ^ 2 ) = 0, with 

(5.2) w{ti,t 2 ) =tiw\+t 2 W 2 + 0{\{ti,t2)\^) in X as {ti,t 2 ) ^ t). 

In a neighborhood of (0,0*, A*) in O x (—oo,0) x (0,7r), the union 
S^* U captures all nontrivial solutions of (5.1). 

(a) Let 0 < 5 < 1. // ni I n 2 , there exists a smooth family of nontrivial 
small-amplitude solutions 

5))* := {(w{r,v),a{r,v), X{r,v)) : 0 < r < e, |sin(n)| > (5} 

of (5.1) in O X (—00, 0) X (0, tt), passing through (0, a*, A*) when r = 0, 
with 

w{r, v) = r cos{v)wl -\- r sm{v)w 2 + O(r^) in X as r —)■ 0. 

In a neighborhood of {0, a*, X*) m Ox (— 00, 0) x (0, tt), the set5^*U5^*^^ 
contains all nontrivial solutions of (5.1) such that |sin(u)| > 6 in their 
projection rcos{v)wl -|-r sin(n)r (;2 on ker Du,iF{0, A*) along Xq. 

Proof. The Lyapunov-Schmidt reduction in Lemma 3.13 tells us that there 
is a neighborhood of (0, A*) m. X xlA \n which (2.10) is equivalent to the 
equation <L(ti,t 2 , A) = 0, where 

(5.3) <h(ti, t2, A) := + t2W2 + tfitiwl -\- t2W2 ,A),A). 

Recall that Z = span {«;}, }, where wj := If we let Uj denote the 

projection onto the span of Wj along the image of £(A*), then Hz = Hi -|-n 2 . 



TRAVELING GRAVITY WATER WAVES WITH GRITIGAL LAYERS 


25 


Defining ‘Ilj := for j = 1,2, the equation = 0 can then be 

rewritten as the system of equations 

(5.4) ^ 2 , A) = 0 , 4>2(ti, ^ 2 , A) = 0 . 

Let us first consider case (i), where ni f n 2 - We claim that 

(5.5) d>i( 0 ,t 2 ,A) = 0 , ^ 2 (ti, 0 ,A) = 0 , 

for all ti, t 2 and A. We will show the first identity; the proof of the second 
being similar. Using that ker A*)|j^(n 2 ) = spanjtc^}, an application 

of the Lyapunov-Schmidt reduction in yields a function mapping 

into satisfying 

(5.6) {I - U2)X{t2W*^ + ^{t2wl,k),K) = 0, 

for all t 2 in a neighborhood of 0 G M. Due to the 27r/(n2N)-periodicity of 
J- {tw 2 +^{tw 2 , A.), A), (5.6) holds with 112 replaced by 11^, whence uniqueness 
of ijj yields 'ip{t 2 W 2 ) = ip{t 2 W 2 )- The hrst identity in Equation (5.5) then 
follows by definition of <I>i. 

Following the proof of Theorem 4.3, we now introduce the function Ti by 

(5.7) Ti(ti,t 2 ,A) ;= / ^iti{zti,t2,A)dz, 

Jo 

and similarly the function 'I' 2 . Then = tj^j by (5.5), and we get different 
cases for (5.4) depending on the values of ti and t 2 - If both ti and t 2 vanish, 
we get trivial solutions. When t 2 = 0 but ti ^ 0, we know from (5.5) that 
the system reduces to <hi(ti, 0 . A) = 0, the solutions of which correspond to 
5 (^ 1 )^ while when U = 0 but ^2 7 ^ 0 we get solutions of ( 2 . 10 ) that he in 
5 (^ 2 ), xhe remaining case is where both ti and t 2 are allowed to be nonzero, 
which amounts to investigating the solutions of 

'I'i(ti, t 2 , A) = 0 , 'L 2 (U, ^ 2 , A) = 0 

in a neighborhood of (0,0, A*). We will look at those solutions with n = . 

For this, we intend to use the implicit function theorem. Note that from 
(5.7) we get ^^(OjOjA*) = (0,0, A*), and that by definition (5.3) of 

we have 


{tut2,A*) = UzD^Xiw + A*),A*){w* + A*)w*), 

where w := tiwl + t 2 W 2 - Using the properties of ip and Theorem 3.11, we 
can therefore conclude that Tj(0,0,A*) = 0. It remains to show that the 
derivative ;^)('Fi,'F 2 )( 0 ,0, A*) is invertible. 

Using analogous computations to those after Equation (4.7) in the proof 
of Theorem 4.3, we have 


^jfi{0,0,A*) = U,DpC{A*)(l)j 


{DpC{A*)4>^,w*)y 


-w. 


w 


’j\\Y 


for j = 1, 2 and /3 = a, X. It follows that ;^)('I'i, 'I'2)(0,0, A*) is invertible 
if and only if the determinant 


{DxC{A*)cj)l,wl)y {DxC{A*)cj)*2,W*2)y 
{D^C{A*)<Pl,wl)y {D^CiA*)<p*2,w*2)y 
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is nonzero. The inner products appearing in this determinant have already 
been computed in Lemma 4.1 and Lemma 5 . 1 . Using these, and elementary 
properties of determinants, we have 

2 / \ \ 2 


C = A 


sinh(6'„i' 


0. 


ni 




sinh(0, 


n2, 


^712 


f{0: 


ni. 


Observe that of 9 ni and 6n2, only 0^2 can vanish. Suppose for the moment 
that also 0^2 7^ 0 . Then 




TT 




TT 


sinh(0„ 


sinh(6*,i 


duj 


sinh(0„ 


— Or, 


cosh(0„/ 


sinh(0„ 




-3 '■ 'O 

whence the determinant in Equation ( 5 . 8 ) can be written as 

2 / / /I \ \ 2 


^ _ Try! / sinh(6*nj ^ ^ 
2 k V On^ 


■'n2 


r(A-)(r(A-)-l)(T_ 1 I, 


^n2 


02 


where A is nonzero due to the assumption of transversality. Hence, we 
immediately see that C is nonzero if and only if r(A*) ^ { 0 , 1 }. A similar 
computation shows that 




6k 


0 , 


ni 


/O 


when 0^2 = 0 . This concludes the proof of part [i). 

Next, we move on to case (ii)^ where ni | n2- We still find that 

^i(0,t2,A) = 0 

for all t2 and A, and so we can introduce 

'I'i(r, u, A) := / (zr cos(u), r sin(u). A) dz 

Jo 

as before; only now written using the polar coordinates (iden¬ 
tifying C and M2^_ Then d>i = For $2, the corresponding identity 

in ( 5 . 5 ) is no longer true in general, but we still have <1>2(0,0, A) = 0 . We 
therefore introduce ^2 through 

T2(r,-u,A) := f [^2t3{zre'''^,A)cos{v) + ^2t2{zre'''",A)s\n{v)]dz, 

Jo 

which yields <1>2 = r'^2- 

Like for case (i), the solutions of d>(0,t2,A) = 0 near (0, A*) for t2 7^ 0 
correspond to solutions in When ti ^ 0 , also r 7^= 0 , and so ( 5 . 4 ) is 

equivalent to the problem 

Ti(r,u,A) = 0, 4'2(r,u, A) = 0, 

which we will now consider. Again, we will use the implicit function theorem 
to find solutions with fj, = ^*. Due to similar computations as those for case 
(i), we have 'l'i(0,u, A*) = 'I'2(0,u, A*) = 0 and 

.T, m A*A {Dl 3 CiA*)cPl,wl)y 
Tl^(0,U,A ) = - 


\w 


Illy 


( 5 . 9 ) 



TRAVELING GRAVITY WATER WAVES WITH GRITIGAL LAYERS 


27 


for all V and (5 = a, \. To find the derivatives of 'I' 2 , note that 

^ 2 / 3 ( 0 , V, A*) = $ 2 ti/ 3 ( 0 ,0, A*) cos(i;) + d> 2 t 2 / 3 ( 0 , 0, A*) sin(r;), 

and so 


(5.10) 


w. 


2iiy 




m 


2llY 


where we have used that cos(niKx) and cos(n 2 Nx) are orthogonal in L^(M). 

From the preceding, we see that the derivative Dj-q, 'I' 2 )( 0 , "u, A*) is 

invertible if and only if the determinant 


(5.11) 


C := 


{DxC{A*)ri,w\)y sin{v) 

{Do,C{h*)(t)\,wl)y {DaC{A*)(j)*2, wl)y sin(i;) 

= C sin(r;) 


is nonzero, where C is the determinant introduced in (5.8). We know that 
C 7 ^ 0 under the assumptions of the theorem, so we can apply the implicit 
function theorem at (0, v) if sin(u) 7 ^ 0. This can be done uniformly in v as 
long as sin(r;) is bounded away from zero. □ 

Remark 5.3. In case (i), the surface profiles in 5^* \ have multiple 

crests and troughs in each minimal period, at least when (ti, t 2 ) is sufficiently 
small. This follows from the asymptotic formula in (5.2). 


Remark 5.4. Observe that the second special case listed in Example 3.7 has 
ni I n 2 , while the third has ni f n 2 . They therefore fall into different cases in 
Theorem 5.2. 


5.1. Properties of the bifurcation sheet. We will now present some 
properties of the sheets of solutions that were found in the two-dimensional 
bifurcation result. Theorem 5.2, following the lines of Section 4.1. The main 
purpose of this is to show that these sheets, found by bifurcating with respect 
to A and a, do not, in general, coincide with the sheets found in [11, Theorem 
4.8]. Like for one-dimensional bifurcation. Theorem 5.2 differs from the one 
in [ 11 ] by the use of A instead of and the addition of the transversality 
condition (4.2). 

The first step towards showing that the sheets differ is Proposition 5.5, 
which is the two-dimensional counterpart of Proposition 4.4. 

Proposition 5.5 (Gradients of a and A). For the solution sheets obtained 
in Theorem 5.2, we have the following: 

• In case (i), the solutions satisfy 

Va(0,0) = VA(0,0) = 0. 

• In case (ii), we have 

ar{0,v) = Ar(0, u) = 0, 

as long as n 2 7 ^ 2 ni. 
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Proof. The proof for case (i) is a simpler variant of that for case (ii), so we 
focus on the latter. By definition of a and A, we have the identity 

'^j{r,v,n*,a{r,v),X{r,v)) = 0 

for j = 1,2. Through taking derivatives with respect to r, this implies that 

(5.12) ^jr + ^jaar{0, v) + TjAAr(0, v) = 0, 

where the derivatives of Tj are evaluated at (0, r;, A*). This linear system of 
equations can be solved for ar{0,v) and Ar(0, r;) because the determinant C 
in Equation (5.11) is nonzero. In order to show that the derivatives vanish, 
it is therefore sufficient (and necessary) to show that '£'^>(0, u. A*) = 0 for 

i = i,2. 

Using the definitions of Ti and 'I' 2 , we find 
Ti^(0,i;,A*) = ^cos(u)d>itjij(0,0,A*)+ sin(u)<l>it^t2(0,0,A*), 
T2r(0,u,A*) = ^cos2(u)$2titi(0,0,A*) + ^sin2(u)d>2t2t2(0,0,A*) 

+ sin(u) cos(u)<h2qt2(0, 0, A*), 

where 

(5.13) {DlT{0,A*){w*,w*),w^)y 

= - 1|',||2 - 

HWy 

for i,j,l = 1, 2. Using orthogonality in L^(M), like in the proof Proposition 4.4, 
one can show that the derivatives in Equation (5.13) are zero, except possibly 
when ni = rii + nj or = In* — nj\. This is only the case when n 2 = 2ni 
and either i = j = 1 and I = 2 or i ^ j and 1 = 1. □ 

We now show that Q'r.(0,u) and Xr{0,v) can indeed be nonzero when 
n 2 = 2ni, which is not covered by Proposition 5.5, by considering the second 
special case listed in Example 3.7. 


Proposition 5.6 (Special case). Let a be the smallest positive solution of 
xcot(x) = 1. When n = (t/\/3 and A* = (1/(2k), — 4K^,7r/2), we have 


^lr{ 0 ,V,A*) 


^2r{0,V,A*) 


1 + ) sin(ti)t()i, 

O 

— + - cos^(cj) ) cos^(u)rt; 
16 2 ^ ^ ^ ^ 


* 
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and so by (5.12), (5.9) and (5.10) that 


ar(0,u) 

= M 

Ar(0,i;) 



(1 + 1 ^ 2 ^ 


sin(u) 

\ cos^(it)) cos(u) cot(u)J 


for a nonsingular matrix M not depending on v. In particular, 0^(0, u) and 
Ar(0,u) are both nonzero, except possibly for isolated values ofv. 


Proposition 5.6 shows that the sheets obtained in Theorem 5.2 are, in 
general, not the same as those obtained in [11]—at least when n 2 = 2ni. 
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5.2. Local description of solutions. We finish by using [11, Theorem 4.8] 
to prove a two-dimensional version of Theorem 4.10, describing all nontrivial 
solutions in a neighborhood of a point falling into case (i). Let therefore 
A* G ^ be such that M(A*) = {ni,n 2 } with ni < n 2 and ni f n 2 , and such 
that either r(A*) ^ {0,1} or 6 {n 2 , a*) = 0. 

Proceeding as in Section 4.2, we use Theorem 3.4 to conclude that there is 
a neighborhood of A in which 


M(A) 


{ni,n2} a = a*,/U =/U*(A), 

< {nj a^a*,^ = ^i{a,X), 

0 otherwise. 

\ 


We may now apply [11, Theorem 4.8] to each point on the graph of /U* (where 
r(A) = r{A*) and a = a*), obtaining a family S of bifurcating solution sheets. 
In addition, one has the solutions in and which were described 

before Theorem 5.2. These are all the nontrivial solutions near (0, A*). We 
omit the proof, which is essentially the same as for Theorem 4.10. 


Theorem 5.7 (Local description). The family S of solution sheets bifurcating 
from points (0,//*(A), a*, A) for X in a neighborhood of X*, together with the 
families and constitutes all nontrivial solutions in a neighborhood 

of{0,A*) inXxU. 

Remark 5.8. We could alternatively have used Theorem 5.2 at points where 
the transversality condition is fulfilled. It follows that, locally, the same 
solutions can be found through bifurcation with either n or A. 


Appendix A. Derivatives of W 

The purpose of this appendix is simply to record the derivatives of J- 
with respect to w at (0, A), up to the third order. These are used to obtain 
derivatives of the bifurcation curves from Theorem 4.3 and the bifurcation 
sheets from Theorem 5.2. 

We have 

D^,T'i{0,A)w = (1 - 'ipQs)'!] + 'ipQs^s 

K)w‘^ = 3V’Lr?^ + 'iplsVl - ^f’OsTl^s + 

DIJ='i{0,A)w^ = -12V’os??^ - 

for J^i, where we suppress the evaluation at s = 1, and the derivatives 
D^X 2 { 0 , A)rc = -2ipossr] - stposTj^x + (d^ + - a)(j) 

D^J'2(0, A)u;^ = + 2sil;osr]r]xx + (dsiAos + ‘^s‘^i^0ss)r]l 

AtJ^ss ^STjx^xs 

dI^J='2{0, A)w^ = - 24ilj0ss'n^ - Gsifosv’^Vxx - {24:Sifos + 12s^V’0ss)Wx 

18r]‘^$ss + 1-2sr]r]x^xs + Qsrir]xx(f>s + 12sry^(^s 

6s‘^rjlfss 


for I' 2 - 
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